The calculation of parameters for Stevens Hamiltonian in a crystalline field generated 
by the electric charge uniformly extended in one direction 



O 

o 
o 

o 

Q 

(N 

O 
O 

Oh' 



I 

C 

O 

o 



> 
o 
o 

(N 

o 
o 

-I— > 
c3 



o 
o 



X 



V. BobrovskiiJ^I A. Mirmelstein, A. Podlesnyak, and I. Zhdakhin 

Institute for Metal Physics, 
Ural Division of the Russian Academy of Sciences 
620219 Ekaterinburg, Russia 
(Dated: February 1, 2008) 

Exact expressions for the parameters of Stevens Hamiltonian are derived within the framework of 
a specific model that assumes uniform character of charge density distribution in a certain direction 
over crystalline lattice. 

PACS numbers: 



FORMULATION OF THE MODEL 



The well-known point charge model is widely used both 
for description of crystal electric field (CEF) effects and 
for analysis of the related experimental results. A stan- 
dard form of its application is based on the so-called 
equivalent Stevens Hamiltonian (ESH) @, |, g, |]. The 
respective coefficients for the Hamiltonian have been tab- 
ulated already for many crystal structures. It should be 
noted that this Hamiltonian is not directly relevant to 
the point character of the electric charges forming CEF. 
Therefore, methods to calculate the Hamiltonian coeffi- 
cients can be generalized for a more complicated case. 
For instance, in 0] they are determined with taking into 
account of screening effects by using of ligand electro- 
static potential in the form of Yukawa. The numerical 
methods for calculation of the ESH parameters are also in 
use. In particular, it is possible to advance point charge 
model for more realistic assuming about electron density 
distribution by introducing a net of negatively charged 
points in interion space ||. (This approach will be de- 
scribed in details in item 6 of the present paper.) 

The perovskite-related compounds comprising Cu02- 
planes in the structure were a subject of intcncc studies 
during the last decade. It is of a general believement 
that the charge state of these planes is responsible for 
the unusual properties of the layered copper perovskites 
as, for instance, the high-temperature superconductivity 
is. The inelastic neutron scattering has proved to be a 
powerful tool for studying the charge states in the high- 
temperature superconductors containing rare-earth ele- 
ments with f-shells {j| ||, ||] . A numerical analysis of the 
data collected in the scattering experiments enabled us to 
obtain earlier some evidences in proof of extended, stripe- 
like charge structures in the Cu02-planes . When an- 
alyzing these data we have found rather general method 
to obtain exact solutions for ESH parameters of a crystal 
field generated by the charge system with density being 
uniform in one direction. 



The specific goal of this paper is the systematic de- 
scription of the mathematical apparatus developed for 
analytical calculation of ESH parameters for this geom- 
etry of charge distribution and also the analysis of the 
consequences resulted from the model. 

Since the charge systems described by our theory are 
assumed to have sufficiently high symmetry (at least or- 
thorhombic) it is convenient to use the local coordinate 
system (x, y, z) shown at Fig. [I] with the origin at the 
rare earth ion and the electric charges situated within 
two parallel planes z — ±H. The density of the uni- 
formly extended in i direction charge can be described by 
the function f{L) where L is the distance in y-direction 
between z-axis and the partial charge filament |l5| ] . 

At first let us consider the electrostatic potential gen- 
erated by the infinite uniformly charged filament oriented 
along x-axis. If the rare-earth ion has f-electron at x, y, z 
then the perturbing crystalline potential generated by the 
filament will be: 



U{x,y,z) = U{y,z) = £ln- 



H 2 + L 2 



{L-yf+{H-~z) 



z\2 ' 



(1.1) 



where £ is the linear charge density. 

Taking into account that due to the small radii of the 
f-shell the inequalities \y\;\z\ <C \/L 2 + H 2 are satisfied 
and introducing the notations: 



V 



M = —- 



y/L 2 + H 2 
2L 



y/L 2 + if 2 ' 
2H 



VL 2 + H 2 



N = — 



VL 2 + H 2 ' 



(1.2) 



one can present the expression (Id) as a polynomial ex- 
pansion: 



(1.3) 



As a consequence of the selection rules for the ESH ma- 
trix elements [Dj, the non-zero contribution into ESH is 
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produced only by the members of expansion (1.3) with 
p + q = 2, p + q = 4 and p + q = 6, for which after 



2 



calculations one obtains: 



where 



F20 — 

Fq2 = 

Fu = 
F40 = 

Fq4 = 

F31 
F%3 

F22 
Feo 

Foe 

F 5 i 
Fib 

F33 
F42 
F24 



2MN - M 3 N, 
2MN-MN 3 , 

-1 + M 2 + N 2 - -M 2 N 2 , 

3 2 6 

i - -iV 2 + N 4 - -N 6 , 
3 2 6 ' 

-3MN + 4M 3 N - M 5 N, 
-3MN + 4MN 3 - MN 5 , 
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-6MN + AM N + AM*N - — M a N a , 

o 

3 5 
1 - -N 2 - 3M 2 + M 4 + 6M 2 N 2 - -M 4 N 2 , 

3 5 
1 - -A/ 2 - 37V 2 + N 4 + 6M 2 N 2 - -M 2 N 4 . 

(1.4) 



1 

- 1" 2 - 




U 6 (y,5) = 










92(L,H) - 


MN, 








- y. 


9i(L,H) = 






9a{L,H) = 



2tig 2 {L,H){~z 2 -y 2 ), 
2£g 4 (L,H)(y i -6y 2 ~z 2 + ~z i ), 
2t;g 6 (L,H){Cz e ~y 6 ) 
-15y 2 ~z 2 (~z 2 -y 2 )}, 

H 2 -L 2 
(i? 2 + L 2 ) 2 ' 
1 H A - 6H 2 L 2 + L 4 
2 



(H 2 + L 2 ) 4 ' 
1 (H 6 - L 6 )-15H 2 L 2 {H 2 
3 



(1.6) 
(1.7) 

(1.8) 
(1.9) 

(1.10) 



L 2 ) 



( J ff 2 + i 2 ) 6 

(1.11) 

The coefficient 2 is introduced into expressions ([T(])-([l^) 
to emphasize that the field in question is assumed to be 
generated by the pair of the charged filaments. Thus, the 
factors (72, 9i, 96 describe the contribution from the one 
charged filament in given symmetric configuration. For 
the system of the charged filaments with the coordinates 
(L,H) (—L,H), {L,-H), (-L,-H) the coefficient 2 in 



(L6)-(L8) is re place d by t he factor of 4. 

The results ( |l.6[ )-( p~TTl ) can be generalized for the 
charge distributed in the plane z = H with the sur- 
face density f(L) which is a constan t in x -dir ection and 
f{L) = f{-L). Then in the formula (p^)-(|T?7|) the linear 
charge density £ will be replaced by f(L)dL. As a result, 
the final expressions for the summary contribution from 
this charged plane in the crystalline field potential in the 
point (x,y,z) look as: 



U 2 (y,S) = 2G 2 -(S 2 -y 2 ), 



(1.12) 



U 4 (y,z) = 2G 4 -(y 4 



6y 2 z 2 



(1.13) 



To obtain the potential generated by the whole system 
of the equally orient ed fi laments it is necessary to sum 
partial contributions (1.3) with their corresponding coor- 
dinates L, H . Below we restrict our consideration by the 
charge systems constructed on the basis of the pair of 
parallel filaments symmetrically situated relatively the 
plane y — (Fig. [|). In particular, it includes such 
important cases as the crystals with orthorhombic and 
tetragonal symmetry. In the local coordinate system pre- 
sented at Fig. [I] our pair has the coordinates as (L, H) 
and (—L,H). It i s obvious that the summation of the 
contributions ( |l.3|) for this pair results in vanishing of 
the terms in which M occurs to an odd power. As a re- 
sult, the relevant part of the electric potential generated 
by the pair of the uniformly charged parallel filaments 
will be: 



U(y, 



U2{y, z) + Ui{y, z) + U 6 (y, z), 



(1.5) 



U 6 (y, z) = 2G 6 ■ {(i 6 - f) - 15fS 2 (S 2 - y 2 )}, (1.14) 



where 



G n =G n (H)= dLf(L)-g n {L,H). 



(1.15) 



In the case of two parallel identically charged plan es sit - 
uated at z = H and z = —H, the coefficient 2 in ( 1.12 )- 
(1.14) is replaced by 4. The expressions (1.12)-(1.14) 
stay formally valid not only for two dimension charge 
but for the volume one too if its distribution described 
by the density f(L,H), i.e. if it does not depend of the 
x -coordinate and f(L,H) = f(—L,H). For instance, it 
takes place when the electric field is generated by two 
cords of arbitrary cross section form. For such a case: 



G n = J dL J dHf(L, H) ■ g n (L, H) 



(1.16) 



and again if the charge distribution is symmetric rela- 
tively not only y = plane but also z — plane, the 



3 



coefficient 2 in (1.1 2)- (1.14) should be substituted by 4 
and 



G n = dL 
Jo 



dHf(L, H)g n (L, H). 



(1.17) 



As in the case of the point-charge model the determi- 
nation of the perturbing Hamiltonian is the evaluation 
of the appropriate electrostatic potential. One-particle 
Hamiltonian for the f-clcctron in a crystalline electric 
field: 



h = eU(y,z). 



(1.18) 



The calculation of the ESH parameters becomes much 
more convenient when using the spherical coordinates 
r, ip in expression for U (y, z) instead of the rectangular 
coordinates x, y, z: 

x = r sin i? cos y = rsini9sin(^, z = rcosi9. 

(1.19) 



On proceeding of this substitution in (1.5), ( 1 .12 )-( 1. 14 ) 
one can obtain after algebraic transformations the follow- 
ing expression for the electrostatic potential: 



(1.20) 



where n = 2, 4, 6; m = 0, ±2, ±4, ... ± n and the coeffi- 
cients A™ are equal: 



if 



4° 



2G 2 -2 



= 2G 2 - 



1/2 



1/2 



2G 4 -±( 7 r) 1 /2, 



.4 



±2 



.4 



±4 



2G 4 



2G 4 



2 /2tt 



1 /2tt 



3 V 35 



1/2 



1/2 



(±2 



1/2 



2G5 . 15 (^)" 2 . 
6 \2730J 



A 



±4 



.4 



±6 



2G 6 



2G 6 



1 /14tt 



7 ( 1 : 



1/2 



1 



231 V 13 



2317T 



1/2 



(1.21) 



It should be noted that during these manipulations we 
used the notations for the spherical harmonics from |Q. 

In real crystalline structures the potential of electric 
field are formed simultaneously by the ligand charges (de- 
scribed, for instance, in framework of the point-charge 



TABLE I: The set of the coefficients A™ corresponding to the 
charge distribution configurations shown at the Fig. ti. 



(a) 



A° 2 4G 2 2(f 



*2 

Ar 

Al 

Ar 

K 

A±2 

A± 6 



1/2 

L/2 



(b) 



4G 2 (2f 
4G4IW 17 
4G 4 |(^) 1/2 

,1/2 



1/2 



4G4 I (If 

4G 6 2(^ 
4G6l5(^) 1/2 
4G 6 i(^ 1 



4G 6 ^r 



2iL 



l 23171 
_L3_ 



1/2 



4G 2 2(f) T7!r 




4G 4 (If 

4G 6 2(£ 


4G 6 i(^; 





(c) 



8G 2 2(f 


8G 







1/2 
1/2 

1/2 



8G4 3 I 35 

8G 6 2 (^ 


Q/"» 1 / 147T 

8G 6 y ( — 




1/2 

j 1/2 
,1/2 



,V2 



model) and by the e xtend e d cha rge structures described 
by the expressions ( 1.2C )-( 1 .21 ). The calculation of a 
such superposition potential and analysis of summary 
system properties can be carried out with the most con- 
venience using the coordinate system associated with the 
crystallographic symmetry axes. Let us consider the case 
when the local coordinate system x, y, z is oriented rela- 
tively the crystallographic system x, y, z as at the Fig. [|. 
It is the orientation that is of the most importance for 
description of orthorhombic and tetragonal systems. In 
the spherical coordinate system connected with the crys- 
tallographic axes: 



x = r sin $ cos ip, y = r sin 1? sin <p, z = rcosi9. 

(1.22) 

Corre spondingly in the local coordinate system (see 
(|l~T9|))we have: 



§ = (p = ip- A. 



(1.23) 



As a result, in (1.20) the replacement arises: 



2m Am 2m —imA 

A n — * A n — A n e 



Then 



U{r,4,<p) = *£A™r n Y™(#,<p) 



= ^i™-e- imA -r"- ^(0,^,(1.24) 

n,m 

where the potential U(r, <p) is assumed to be notated in 
the crystallographic coordinates whereas the coefficients 
A™ are calculated in the local one. Three charge configu- 
rations with different symmetry are shown at the Fig. |^. 
Corresponding set of the ^^-coefficients are presented 
in the Table [| It should be noted that at A ^ 0, § the 
coefficients A™ have imaginary parts. 



2. CALCULATION OF THE ESH PARAMETERS 

Following the procedure proposed in |^] for determi- 
nation of the ESH parameters we shall turn to tesseral 



4 



harmonics from spherical ones in fll.24| ). Accordingly || 
for even m the tesseral harmonics can be expressed via 
spherical harmonics as: 

Z n = Y n i 



Zi 



f ym I y— ml 
^2 in n J ' 

i^ 1 n " 



(2.1) 



where m=2,4,6 m < n. Accordingly: 
1 



V2 
1 

V2 



S 7 C — \7 S X 
X nrn nmJ' 



(2.2) 



Takin g into account (2.1)-(2.2), one can obtain from 
(pi: 



(2.3) 



where a = 0, c, s. 

Then for the charge system formed on the basis of the 
set of parallel filaments (i.e. with the potential de scrib ed 
in the local coordinate system by the expressions (1.21)), 
we have in the crystallographic coordinate system: 



7° 

7nm 
Tnm 



V2i"cos(mA)-(l-<5 m0 ), 
V / 2i™sin(mA) • (1 - S m o). 



(2.4) 



There exists the universally adopted method for getting 
rid of s-harmonics Z^ m by the corresponding choice of 
the coordinate system where 7* TO = 0. However, aiming 
to derive formulas suitable for calculation of superposi- 
tion potentials we shall not restrict extension of our ex- 
pressions. The collection of formulas ji] expressing the 
tesseral harmonics as a functions of the rectangular coor- 
dinates is completed below with the s-harmonics. Using 
the notations: 



K„ 



K„ 



■I 1 1 



Z' = 



fO 



(2.5) 



and the expressions ( 1.22 ), ( |2.l| ) and also the spherical 
harmonics from |lj], we have obtained the formulas for 
the f„ m ; /* m listed in the Table ||. 

The Hamiltonian of the f-shell in the CEF is as follows: 



H c f = 



(2.6) 



where we are normally concerned with the summation 
over the f-electrons. Thus: 



Hd = -|e| £ £ £ 7nm ' K nm ■ £ W, Zi), 

ri m > ol i 

(2.7) 

Using Wigner-Eckart theorem (l], |2) one can replace: 



Then: 



where 0° 



(2.8) 



(2.9) 



n m>0 a 

operator equivalents introduced by Stevens 



and 



B nm — \e\ r y nm K nm Q n {r ), 



(2.10) 



The so-called irreducible matrix elements 0„ are listed 
in (!]]. They are often denoted as 9„ = aj;[3j;^j for 
n=2, 4, 6, respectively. The averaged value for the nth- 
power of the f-she ll radii is denoted as (r n ) Taking 
into account (2.4), ( 2.10| ) let us write down now the ex- 
pressions for the B" m for the charge system characterized 
by the angle of rotation A relatively crystallographic co- 
ordinate system as at Fig. [|: 



B 20 

B 22 
B 22 

B 40 
B 42 

B A2 
B ti 
B li 
B 6Q 



B 62 



\b2, 

icos(2A)6 2 , 
isin(2A)& 2 , 



-cos(2A)6 4 , 
-sm(2A)& 4 , 
- cos (4 A) b 4 , 

o 

1 



sin(4A)6 4 , 



= u 



B 62 = 



B 64 - 



B 



64 



B, 



(id 



B, 



16 
15 

32 
15 
32 
_3_ 

16 

3 

ic 
1 

32 
1 

32 



cos(2A)6 6 , 
sin(2A)6 6 , 
cos(4A)& 6 , 
sin(4A)6 6 , 
cos(6A)fo 6 , 
sin(6A)6 6 , 



(2.11) 
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TABLE II: Components of the tesseral harmomics expressed in Cartesian coordinates. 
K f c f 

J-^jim J n?n J nm 



2,2 Hv) 1/2 x2 -y 2 2x v 

4,0 ±{\) 1/2 35z 4 - 30z V + 3r 4 

4,2 1/2 (7z 2 -r 2 )(z 2 -y 2 ) (7. 2 - r 2 )2zy 

4,4 ^(f) 1/2 ^-G^V+y 4 4^ 2 -y 2 ) 

6,0 i (13 ) 23l2 6 _ 315 ^4 r 2 + 1Q5z 2 r 4 _ 5r 6 

6,2 i (4 s ) V2 {16^ - 16(x 2 + y V + (x 2 + y 2 ) 2 }^ 2 - y 2 ); {16z 4 - 16(^ 2 + y V + (z 2 + y 2 ) 2 }2xy 

M §(#) 1/2 (Hs 2 -r 2 )(s 4 -6xV + 2/ 4 ) (llz 2 - r 2 )4*y(* 2 - y 2 ) 

6,6 fH^) 1/2 ^ ~ 15xV + lSzV - y 6 2xy{4( a: 2 -y 2 ) 2 -(x 2 + y 2 ) 2 } 



2,0 if!) 172 3 Z 2 



where: 



Correspondingly: 



b n = -C\e\G n e n (r n ), 



(2-12) 



and the coordination number C is equal 2 for the charge 
distributed within the z = H plane with the density 
f(L) = /(— L), and C = 4 when the charge is placed 
within the z — ±H planes with the density f(L,H) = 
f(-L,H) = f(L,-H) = f(-L,-H). Naturally, the co- 
efficient C takes the s a me va lues for the cases described 
by the formulas ( |1.15| ), (1.16) respectively. 

The Stevens operator equivalents are listed for a — 0; c 
in [jjj. The expressions for the operators with a = s can 
be derived by means of well-known procedure consist- 
ing of the replacement of the coordinates x,y,z in the 
expressions for /" m by the corresponding components 
J x ',Jy',Jz of the angular momentum operator J. It is 
the reason why the formulas for the tesseral harmonics 
in CEF-formalism are transformed back to Cartesian co- 
ordinates. It should be noted that it is conventional in 
Stevens formalism to use the operators J a without h in 
the commutation rules, i.e.: 



J. x , J y i J z 



Jy\ Jz 



J zi Jx — i-J y , 

For the operators J± defined as: 

J± = Jx ± iJy, 

we have: 



(2-13) 



(2.14) 



J±\J, M) = y/JiJ + 1) - M(M ± I) • | J, Mil), 

(2-15) 

However, as long as the commutation rules for the op- 
erators x;y; z and J x ] J y ] Jz are different, the replacement 
procedure demands preliminary symmetrization of f^ m 
over all possible permutations of the coordinates x, y, z. 
For /| 2 , /| 4 , /| 6 this manipulation is trivial. For example: 

& = 2x y = ji {( x + w) 2 - < x - '"I 2 ) - ( 2 - 16 > 



And: 



/44^{(z + iy) 4 -(*-h/) 4 }, 



°- = i{^- J -}' 



/66 = ^{(z + iy) 6 -(*-i2/) 6 }, 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



a 



66 



^{J 6 + -J 6 -}, (2.21) 



However, the method that requires less tedious alge- 
braic manipulations than the universally adopted one, 
may be applied to obtain the expressions for arbitrary 
operators 0" m . Let us take into account that: 



f s 

J nm 



l l 

Knm v2i 



7 (r n Y™ - r n Y~ m ) 



(2.22) 



The functions Y™; Y^~ x \ ■■■Y r T n themselves form the full 
set of the irreducible tensor operators. Thus, accordingly 
Wigner-Eckart theorem the combinations f-Y^ 1 may 
be directly (i.e. without using of rectangular coordinates) 
substituted by the set of the equivalent tensor operators 
of the same rank: 



>?y-(^,^) = e n <r")T^, 



(2.23) 



Accordingly the algorithm proposed in |12[ , one can 
choose the operator proportional to J™ as the equivalent 
tensor operator W£. Other operators of the full set can 
be obtained with using the recurrent formula: 



J_;VF™ = y/n(n + 1) - m(m - TjW™- 1 , (2.24) 
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It is easy to see that by accepting the definition W£ = 
-^K nn J™ we completely match this procedure with the 
formulas ( pTl7| ), ( pTl9| ), so that: 



b( 



\/2,K nm 

i 



T4 7 " 



The calculations for the 0* m yields: 
1 



rt n 

w~ 



(2.25) 



0\ 2 = -{J 2 [7J| + 14J Z - J(J+l) + 9]- 



Ly 62 — 



J 2 [7J Z 2 - 14J Z - J(J+l) + 9]}, (2.26) 



-{J 2 [33J 4 + 132 .11 - 18J(J+ 1)J 2 + 

273J Z 2 - 36J(J + 1) J z + 282J 2 + J 2 (J + l) 2 - 

26 J( J + 1) + 120] - J 2 [33 J 4 - 132 jf - 

18 J( J + 1)J 2 + 273 J 2 + 36 J( J + 1)J Z - 

282 J z + J 2 ( J + l) 2 - 26J( J + 1) + 

120]}, (2.27) 



O, 



04 



{ J 4 [11 J 2 + 44J Z - J(J + 1) + 50] - 
Jl [11J 2 - 44 J z - J{J + 1) + 50] } , (2.28) 



2i 



Jointl y with t he fo rmula s list ed in 1, the expressions 
( ^17| ), ( |2~19| ), dgjgjj ) and (p^-(j2.28|) give the full set of 
Stevens operators for n=2, 4, 6 an d e ven m . It should 
be noted that our expressions ( 2.26 ) - ( |2.28| ) contain the 
difference of two mem bers corresponding to the opera- 
tors W™ and W~ m in ( [2.25 ) respectively. Replacing this 
difference by sum and the factor l/2i by 1/2, one can 
obtain, as it is clear from the algorithm used, the expres- 
sions for the corresponding O c nm . It is easy to make sure 
after simple transformations that these coincide with the 
formulas from Q. 



3. THE DETAILS OF THE CHARGE 
DISTRIBUTION 

Main distinctive feature of the charge systems studied 
in this paper is their extended ( "filament-like" ) charac- 
ter, i.e. the absence of dependence for charge density 
distribution in the x direction. As it can be seen from 
11), this fact has already estated the strong interre- 



(1 



lations between the coefficients B" m at given n. For in- 
stance, B22 = cos(2A)i?2o- 

It should be stressed that these relations are indepen- 
dent of the other details of the charge distribution. The 
similar relations between the coefficients for more com- 
plex superposition systems are demonstrated in the H. 



Let us consider more closely the case when the electric 
charge is distributed within the z = H plane with the 
surface density f(L) = /(— £), i.e. when the coefficients 
G n are determined by ( 1.15 ). In trod u cing th e par ameter 
A = L/H, one can obtain from (|il|-(|l.ll|), fllT5| ): 



G n = C n d\f(HX)xn(X), 
Jo 



where: 



Cn 



1 

X 2 (A) 
X 4 (A) 
Xe(A) 



1-A 2 



C 6 



3H 5 ' 



(3.1) 



(3.2) 



(1 + A 2 ) 2 ' 
1 - 6A 2 + A 4 
(1 + A 2 ) 4 ' 
1 - A 6 - 15A 2 (1 - A 2 ) 
(1 + A 2 ) 6 ' 



(3.3) 



The formulas (3.1)-(3.3)) describe the contributions of 
the "partial filaments" with the coordinates L = HX into 
the potentials of the 2-nd, 4-th and 6-th orders respec- 
tively. The functions x n (A) play the role of the weighting 
factors. It is clear from Fig.|| that these contributions os- 
cillate with A not synchronously. Moreover, the rates of 
their damping are also sufficiently different. 

At last for a certain practical calculations it is rather 
convenient to define the parameter t as follows: 



Then: 



L = fftan(-), (3.4) 



G n = C n / dtf(Httm{-)) ■ fL n (t), (3.5) 



n 2 (*) = cos(i), (3.6) 



fi 4 (i) = - {cos(3i) + 2 cos(2i) + cos(t)} , (3.7) 



fi 6 (i) = — {cos(5t) + 4cos(4t) + 6cos(3t) + 
16 

4cos(2i) + cos(t)}, (3.8) 
The plots for the functions fl n (t) are shown at Fig. ||. 

4. THE MODEL OF THE CHARGED STRIPES 

Let us consider the case when the electric charge is 
distributed along the plane z — H within the system of 
the parallel stripes uniformly charged with the surface 
density Fq. As usual, it is assumed that f(L) = f(—L). 
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The contribution into G n from the stripe situated in the 
range L N < L < (L N + Wm) can be calculated using 
(|3.l[) or (|3.5[). As a result one obtains: 



(4.1) 



^2 (A) 
vi (A) 
fe(A) 

where the variables 



A 



1 + A 2 ' 
A(3 — A 2 ) 
3(1 + A 2 ) 3 ' 
A(5- 10A 2 + A 4 ) 
5(1 + A 2 ) 5 : 



aw 



Pn = 



L 



N 



w. 



N 



H 



(4.2) 



(4.3) 



are the values of the parameter A corresponding to the 
stripe left and right edges. The summary contribution 
from the total system of the charged stripes can be cal- 
culated accordingly: 



(4.4) 



N 



At last it is useful to write down the final formula for the 
coefficients b n taking in to ac count the contribution from 
the N-th stripe. From ( |3| ) and (|yj) it follows: 



= -Ce 2 a(r n )e n C n R(Av) - u n (a N )} , (4.5) 

The symbol a here means the surface concentration of 
the carriers generating the stripe charge, so that: 



Fn 



a e 



(4.6) 



where e is the elementary charge and the variable a is 
positive for the holes and negative for the electrons. It 
is c onv enient during the practical calculations to replace 
in ( flit ) the value e 2 by e 2 = 14400 mcV • 1 A, simul- 
taneously expressing all the lengths in Angstroms . Th e 
coefficient C in (4.5) is of the same sense as in (2.12). 
Carrying out summation of the expression (4.5) over N 
with corresponding values of the angle A in (2.11), one 
can obtain the formulas for the ESH parameters in the 
cases of rather complex superposition systems. 



5. PERIODICAL TWO DIMENSIONAL 
LATTICES OF THE CHARGED STRIPES 

Now we co nside r the case when the surface charge den- 
sity f(L) in ( 1.15 ) is a periodic function with the period 
T. 



f(L + T) = f(L). 



(5.1) 



Introducing the surface density a(L) of the number of 
the carriers: 



f{L) = \e\a{L). 
Expanding u(L) to Fourier series results in: 



. 27rfc 

a{L) = ct + 2^ a k C0S (— L )- 

fc=i 



(5.2) 



(5.3) 



The term <tq in (5.3) corresponds to the uniformly 
charged plane. This part of the charge density generates 
only the uniform electric field with the constant gradient 
and naturally with the zero contributions of the 2-nd, 4- 
th an d 6- t h or ders into the potential. It is easy to see 
from Q-(|3) t hat our re sults satisfy to this condition. 
Then combining (|3.l[) and (5.3): 



G n — \e\C n ^ EnkUk, 

where the coefficients 



(5.4) 



fe=i 



Enk = / d\ x„(A) cos(27rfc— A) (5.5) 
Jo 1 

depend only on the relation H/T and are not connected 
with the other details of the charge di stri bution described 
by the coefficients Uk- The formulas (4J3) stay to be valid 
for the system in question but now they look as: 



b n = -Ce 2 C n e n (r n )J2 E nk<Tk, 



(5.6) 



fe=i 



As it is clear from (\Lq), the full set of the coefficients 
B" m within the framework of the model developed is 
completely defined by the values of the three parameters 
^2, ^4, b 6 , which can be obtained, for example, by fit- 
ting of the experimental and calculated neutron spectra. 
Thus, the experimental information can be a tool for ver- 
ification both of the microscopic and of the phenomeno- 
logical models of the charge distribution with three arbi- 
trary p aram eters. It should be noted that the expressions 
( |5.4|) - (5J3) describe the contributions from one or two 
charged planes dependingly of the coefficient C (C = 2 
or 4, respectively). But they can be easily generalized 
for more complicate configurations. For instanc e, in the 
case of a stack of identically charged planes, in ( |5.6| ) one 
should estate (7 = 4 and substitute the coefficients E n k 
by the E nk : 



E n k — ^ E n k , 



(5.7) 



H 



where the summation of (|j) is carried out over the plane 
coordinates H > 0. 

Now a few words about the periodicity of the charge 
distribution. Well-defined character of the experimental 
neutron spectra gives us the evidence of the fact that 
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the crystal lattice and charge distribution periods are 
commensurate. In other case the spectrum would be 
smoothed. Moreover, the coefficient interrelating these 
periods should be not large due to the same reasons. 
There are three examples shown at the Figj^. 

• The period of the charge distribution is equal to 
the lattice one. In such a case we have identical 
positions for every Rare-Earth ion and only one set 
of the crystalline electric field parameters. 

• If T = 2a, the situation is different. We have two 
types of the REI positions with their own crystal 
field. If the condition f(L) = f(—L) is satisfied 
then both positions are symmetric. 

• If T = 3a there exist three types of positions. But if 
we deal with the systems with at least orthorhom- 
bic symmetry then the charge distribution in the 
elementary cell should be rather symmetric too. In 
other case the distortions seem to be inevitable. 
The condition f(L) = f(—L) guaranties, at least, 
the symmetry of distribution for one type of posi- 
tions. Then immediately we get rather interesting 
thing: there are two different positions ("1" and 
"2") for the REI in such a system. But the "I" po- 
sition is asymmetric. As a consequence there must 
be the B nm coefficients with both even and odd m. 

As a rather general example, let us consider the model 
of two parabolas (Fig||) when the period T = a and the 
function <r(L) in the range of (0, |) is assumed to be 
equal to: 

a(L) =a + kL 2 + Q(L - L )| {s 2 - (L - ^) 2 } , 

(5.8) 

where Lo — | — 5; oo is the constant that does not affect 
t>2, b 4 , b$ and should be obtained by means of a certain 
additional speculations 



e(L-Lo) 



I, L > L , 
0, L<L . 



As a result, the system of the three equations (at n = 
2, 4, 6) in terrelating the variables k, fi, S may be derived 
from (O), 0: 



Pn = KQn + ^S n (S), 



where: 



Pn 



b n TT 2 



Ce 2 ® n (r n )C n a 
(-If 



Qn — a ^ E n k 

k=l 



k 2 



k=l 



(-1)* 

k 2 



2nkS 



sin(- 



(5.10) 
(5.11) 
(5.12) 



,2-Kk8' 

cos( ) 

a 

(5.13) 



The equations (5.10) are linear relatively the variables 



k, fi. The solution of this system easily can be reduced 
to determination of S from the following transcendent 
equation: 

(Q2P4 - Q4P2)S 6 (S) + (Q 6 P 2 - Q 2 P e )S 4 (S) + 
{QiPe-Q & Pi)S 2 {8) = 0. (5.14) 



After solving of ( 5.14 ), other parameters will be deter- 
mined immediately: 



/' 



Q2P4 — Q4P2 

Q2S46) - Q 4 s 2 (sy 

P 2 S 4 (S) - P 4 S 2 (5) 
Q2S i (6)~Q i S2(S)' 



(5.15) 



Sometimes from certain additional speculations (for 
example, concerning the concentration of dopants in 
the crystalline lattice) one can independently indicate 
the concentration of the carriers injected into the plane 
z = H . Let us denote the total number of these carri- 
ers per the element a X d of the charged plane shown at 
t he F ig.|| as A^. In such a case the expressions (5.14) - 
fl5.I5| ) can be completed by the formulas for determina- 
tion of the (Tq term from (5.8). Really: 



Thus, accordingly this model the function cr(L) in the 
range of (0, |) is a superposition of the two parabolas 
with the extremums in the points L — and L = a/2 
respectively. The coefficients k and fi can be of arbitrary 
signs. Also taking into accou nt t he possibility of adding 
of an arbitrary constant to (5.8), one can see that this 



formula describes rather wide set of the charge distribu- 
tions. Aft er c alculation of the Fourier coefficients for the 
function (5.8) we obtain: 



(-If 
cos( 



ir 2 k 2 
2nk5' 



2fj, 



2nkS 



,2Ttk8, 
sm( — )- 



(5.9) 



N Q = 2d dLa{L), 



Substituting fl5.8|) into (5.16) yields 



A 

°o = — r - 77; 



ad 



10 a 2 n- —n8 2 , 
12 6a 



(5.16) 



(5.17) 



Sometimes the calculation of b n is convenient not making 
use of Fourier expansion s. F or example, direct summa- 
tion of the expressions (4^) over N is po ssible for the 
model of the charged stripes. Accordingly (2.12): 



-C-|e|-e r , 



N=l 



(5.18) 
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For the periodical lattice of the stripes charged with the 
surface density F = cr\e\ we have: f(HX) = F at A <E 
(ajv;/3jv) and f(HX) = at A <f_ (cxn;0n), where 



a N = (Na-^^)/H 
Pn = (Na - °-^-)lH, 



(5.19) 



W is the width of a stripe. 

I t is p ossible to accelerate the convergence of the series 
in ( 5. IS ) using the fact that one can add an arbi trar y 
constant to the charge density in the expression (3.1). 
For example: 



G n — C n 



dX- 



f(HX) - 



F W 



Xn(X), (5.20) 



where we have extracted from the charge densi ty its aver- 
age value. In practice, the summation in ( 5.1 8| ) is carried 
out up to a certain maximal stripe number N = M . Rep- 
resenting ( 5.2CD as: 



dX- 



dX- 



f(HX) - 



F W 



f(HX) - 



F W 



a 



■Xn(X) + 

Xn(X), (5.21) 



we then take into account that at large (3m , as it can 
be seen from the Fig|| the second integral contains the 
product of a slow varying function and the function oscil- 
lating around the zero average value. It gives us the pos- 
sibility to neglect the second integral (from the physical 
point of view it corresponds to the mutual compensation 
of the two neighboring stripes carrying the charges of the 
equal value but of the opposite signs). Then calculating 
the first integral we obtain: 



G n w C n Fo 



M 



w 



Eyy 
{v n {M - v n (a N )} v n ((3 M ) 



.N=l 



(5.22) 



Correspondingly: 



b n = -Ce 2 aQ n {r n )C n 



W 



M 



LN=1 



^n(ajv)} - —v n {(3 M ) 



(5.23) 



POINT CHARGE MODEL CALCULATION 
OF THE ESH PARAMETERS FOR THE 
COMPLEX CHARGE CONFIGURATIONS 



Below we list the formulas adopted for numerical cal- 
culation of the coefficients B" m generated by the clouds 
of charged points. Then we use the results obtained for 
evaluation of an applicability of the charged filaments 
model to description of real systems. In the point-charge 
model the coefficients are expressed via the spher- 
ical coordinates Rj , dj , ipj of the charges generating the 



crystalline electric field as M 



E 

3 



2n + 1 



R 



■n+l 



(6.1) 



The necessary expressions for the tesseral harmonics 
are listed in Table [Ey. In this section we use for the 
Cartesian coordinates of the charges the traditional cap- 
ital letters Xj , Yj , Zj instead of Dj , Lj , Hj that are re- 
served only for filaments. As a result, one can obtain 
from ( [2.10 ) and Table || the following list of the formu- 
las: 



B 2 o = -\e\-e 2 (r*)J2lr ' ' 



R 5 - 



B c 22 = -\e\*e 2 (r*)J2<l/> Y > 



R- 



B iQ = -\e\-Q^)Y: 



35Zf - 30ZfRj + 3i?j 
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#42 = -|e|— e^r 4 )^^- i 



42 



e|^e 4 (r 4 )^ gj 



(7Z| - ^ 2 )2X J -y, 



.B44 = 



i i 35 n / *vV^ ^M^J-lll 



1 2312! - 315Z 4 i? 2 + W5Z?Ri - 5i?f 
5(50 = -I e l256 08 < r >E® " ^ — 



^ 2 = -|e|^e 6 (r 6 )5: 



512 



[16Z 4 - ie(x? + y/)z 2 + (x| + y/) 2 ] (xj - y?) 



R 



13 



^ = -|e|— e 6 (r 6 >5> 



[lez 4 - i6(x| + y/)z 2 + (x] + y/) 2 ] 2^ 



Rf 



,,63 n/6 . v (iiz 2 -i? 2 )(x 4 -6x 2 y/ + y/) 

64 = _|e| 256 6 <r ' ^ * 1 15 : 

7 J 



5, 



64 



i|63«/«v\- (H4-fl 2 )4^(X 2 -y/) 

l e l^ 6(^>2^* 



0(5 



I I 231 / 6\ \ " X j _ ^ X t Y J + l ^ X ] Y f ~ *j 

1 '512 bV /Z^^ pis 



B 66 = 



'512 



3 



2x J y J [4(x 2 -y/) 2 -(x 2 + y/) 2 ] 



R 



13 



(6.2) 



The expressions (6.2) are convenient for numerical cal- 
culations in the models representing the electric charge 
distributed in the lattice as a net of the charged points. 
Let us use this formulas for evaluation of the inaccuracy 
induced by the replacement of a fini te f ilament by an in- 
finite one. As it can be seen from (6.2), the coefficients 



decrease with a distance as: 



B 2m ~ R } 



pa p — 5 p 



Dm v j 



(6.3) 



Thus, it is obvious that the replacement in question 
contributes the maximal error at n — 2. So, it is reason- 
able for our evaluation to calculate the coefficients B 2 o 
and -B22 f° r the finite charged filament (exactly speak- 
ing, belonging to the configuration of the filaments in 
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the sense of the formulas ( [2.12 ) oriented along the x-axis 
in the range (-D; +D) and characterized by the coordi- 
nates Y = L, Z = H . Such a configuration automatically 
give s B s nm = 0. Replacing qj — > ^dx and transforming 
(3.2) from the sum to the integral yield: 



Bo 



-M-£-e 2 -(r 2 ) 

i? 2 cos 2 (^ )}, 



1 sin(yo) 

2 ' (H 2 + L 2 )- 



{(H 2 



L 2 )+ 
(6.4) 



= -leK-ea-^.i^.^l^lsin^o) 



L 2 + H 2 



(3-sin 2 (^ ))L 



<y5o = arctan( 



VIP + L 2 



(6.5) 



(6.6) 



For the infinite filament 



Ass 



ummg 



\/H 2 + L 2 <C -D, let us expand the expressions ( 3.4 )-(|675|) 
over this parameter. As a result, the following convenient 
formulas for evaluation of relative inaccuracy appear: 



B 2 a = B 2Q {oo) 



B 22 



1 H 2 



3 H 2 



H 2 



H 2 + L 2 
D 2 



H 2 



2 L 2 -H 2 



D 2 



(6.7) 



(6.8) 



where £?2o(oo)i B^ioo) - are the corresponding expres- 
sions for the infinite filament. Naturally the formulas 
( |6.7[ )-( |6~8[ ) are not applicable for the case L — H when 
-620(00) , i?22(o°) are equal to zero. In the latest case 
only a bso lute evaluations on the basis of the expressions 
(|6.4[)-(|6~5|) are possible. At L = the value of the ratio 



D/H = 10 is follows from (S.7)-( p.8| ) to be enough to 
provide the accuracy of 1% for the model of the infinite 
filament. To achieve the same accuracy at L ~ 5H one 
needs of the ratio D/H ~ 50. However, when the lat- 
tice of the stripes is described, such an increase concerns 
only the relative error of the far stripes input whereas 
the main contribution to the potential is generated by 
the nearby ones. 



7. CONCLUSION 

The new model is developed to strengthen the pos- 
sibilities for analysis of experimental data obtained by 
inelastic neutron scattering. It is particularly valu- 
able for studying the charge distribution in the layered 
perovskite-like compounds doped with the rare-earth ions 
being used as a local probe. In contrast to the widely 
adopted point-charge model, the new approach directly 
takes into account the extended character of the charge 
distribution in one direction. The model provides clear 
analytical expressions for the crystalline electric field pa- 
rameters, which may be used efficiently in further eval- 
uations. As a result, the strong interrelations between 
these parameters were found to originate from the spe- 
cific charge distribution topology. The similar relations 
may easily be established for more complicated, super- 
positional systems. The model enabled us to develop 
methods of the charge profile reconstruction from the ex- 
perimental data. 
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FIG. 1: The local x; y; z and crystallographic x; y; z coordi- 
nate systems for description of the field generated in the point 
r by the charged filaments AB and AlBl. 



FIG. 2: Superpositional charged systems characterized by 
three parameters: (a) - symmetry mmm(D2h); (b) - symme- 
try 42m(Z?2d); (c) - symmetry A/mmm(D ih . 



FIG. 3: Weighting functions x«M f° r the contributions to the 
crystalline field potential of the 2, 4 and 6 orders, respectively. 



FIG. 4: Weighting functions Q n (i) for the contributions to the 
crystalline field potential of the 2, 4 and 6 orders, respectively. 



FIG. 5: Schematic charge distribution profiles with different 
periodicity. 



FIG. 6: The model of two parabolas. The summary carrier 
concentration a{L) consists of the two prabolic contributions 
(7i (L) and <J2{L) and may be added by the constant oo. 
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